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ABSTRACT 


The work presented in this thesis relates to the signal 
processing for NMR imaging The problem of NMR imaging is 
Viewed as a system identification problem The Bloch equation 
governing the NMR imaging is treated as state variable repre- 
sentation of the system to be identified The RF signal exci- 
tation IS treated as the input and the gradient signals are 
viewed as system parameters Thus the NMR imaging system is a 
linear time varying system, and the identification problem 
amounts to finding the initial state given the response (FID 
signal) over an interval of time The existing modalities of 
NMR imaging are explained in this point of view 

NMR chirp imaging, based on stochastic imaging is proposed 
This scheme of imaging exhibits important merits over other 
pulse sequences One is simplicity in generating gradient fun- 
ctions In NMR imaging the quest for high spatial resolution is 
of fundamental concern in clinical applications The proposed 
chirp imaging permits slow switching of field gradients, admi- 
tting stronger gradients for improved resolution In addition 
this scheme may permit, using charge coupled device hardware 
thereby reducing the hardware comolexity 



Computer simulation is carried out for obtaining the FID 
signal for stochastic excitation and chirp imaging The 
appropriate signal processing techniques are simulated and the 
processed images are found to be in close agreement with the 
assumed models 
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CHAPTER 1 


INTRODUCTION 

Nuclear Magnetac Resonance (NMR) is an indispensable tool 
for material characterization in many disciplines Recently 
it has been applied as a method of internal imaging, and in 
this role NMR enjoys advantages like non-invasive , penetrating 
bony structure without any attenuation, non hazardous nature, 
high resolution capability, potential for specific imaging 

In this chapter we discuss the principles of NMR and its 
applications in imaging context The hardware associated with 
NMR imaging system, historical development and system theory 
approach to NMR imaging is presented The NMR imaging problem 
can be viewed as a system identification problem, which is well 
studied in communication and control area and this view point 
IS elaborated in this thesis 

1 1 PRINCIPLES OF NMR IMAGING [l3] 

Nuclear Magnetic Resonance is based on the fact that the 
atomic nuclei of most elements possess a ’spin* which in turn 
IS associated with a small nuclear ’magnetic moment’ The 
nucleus, e g , the hydrogen atom or proton, can therefore be 
imagined as behaving like a tiny bar magnet or compass needle 
(a magnetic dipole) In the absence of an external magnetic 
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field the directions of these dipoles have a random orienta- 
tion (Fig 1 1(a)), but with the application of an external 
static magnetic field the nuclei begin to precess describing 
a conical path around the magnetic field direction with a 
frequency proportional to field strength (Fig 1 1(b)) This 
frequency is called nuclear magnetic resonance frequency or 
Larmor frequency Energy considerations favour orientation 
of the nuclear dipoles in a direction parallel to the applied 
field direction and the majority of dipoles precess such that 
their points tend more towards this direction The sum of 
these elementary dipoles results in a macroscopic magnetiza- 
tion (M) Applying a radio frequency (RF) pulse brings all 
the nuclear spins into phase All the processing dipoles in 
phase, produces an alternating magnetic field capable of 
inducing into a pick-up coil an alternating current of 
resonating frequency which could be measured When the RF 
pulse IS removed, detection of the nuclear signal induced 
into the receiver coil by free precession is then possible 
This induced signal is known as a Free Induction Decay (FID) 
signal 

The resonance phenomenon can be visualized from quantum 
mechanical description of NMR (Fig 1 2) When proton nuclei 
are placed in a magnetic field, the nuclei will be in two 
energy states, namely + jiH and - pH states [?] , where p and H 


(a] 


Fig 1 1 Behaviour ox nuclear spins (a) vvithout magnetic field 
(b) \ ith strong magnetic field 



Fig 1 2 Quantum mechanical description of N^'lR spin 
e xcitation 
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are the nuclear magnetic moment and applied magnetic field, 
respectively For those nuclei in -pH energy state, irradia- 
tion of electromagnetic radiation of energy E equivalent to 
2 pH will raise the proton energy state upto +pH This energy 
IS usually supplied by an RF magnetic field 

If an auxiliary magnetic field is applied by means of 
gradient coils so that, for example, the field of the magnet 
increases linearly from left to right then the protons at the 
extreme left of the body will have the lowest nuclear reso- 
nance frequency, with tho frequency increasing linearly to 
the right depending on the gradient field, thus the spatial 
coordinates are mappod to frequency scale The amplitude of 
the NMR signal at a particular frequency, hence, is propor- 
tional to the number of nuclear spins (eg , protons) on 
resonance at the corresponding spatial coordinate 

To obtain information about the structure of more complex 
three-dimensional form (like human head) having within it an 
uneven distribution of proton concentrations, magnetic field 
gradients need to be applied in all directions so that fre- 
quency values may bo assigned to coordinates throughout the 
object 

The relaxation mechanisms associated with excited nuclear 
spins are transverse relaxation (T2) and longitudinal relaxa- 
tion (Tj^) T2 IS faster than so that T2 is always smaller 
than 
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The characterization of tissue matter is in terms of 
water content Water xs the medium in which most of the body 
chemistry or metabolism takes place, hence it plays an 
essential role in all biological systems [18] The reason 
why water has so far played a minor role in medical diagnosis 
IS simply the fact that until the advent of NMR imaging, 
there really was no way of studying it in vivo Although 
NMR imaging will doubtless find applications in the study of 
plants and fluid transport, emphasis is placed here on the 
study of biological tissue in man Biological tissue con- 
tains on average around 75/ water (Fig 1 3), which can readily 
be detected by NMR methods Quite a large fraction of water 
1 $ contained within the cell cytoplasm but there is also a 
substantial amount of extra cellular water in addition to the 
bulk fluids contained within the body Different healthy body 
tissues contain different amounts of water and variations in 
water content can be used to some extent to differentiate the 
various human tissues although the variations between tissues 
are rather small 

It IS also found experimentally that there are differences 
in the spin relaxation times of various tissues and there is a 
good correlation between the spin-lattice relaxation time 
and water content among the various animal tissues^ the longer 
Tj^'s being observed in tissues with higher water content 
Damadian [8] has shown that T^ is longer in malignant tissue 
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by a factor of about 1 b to 2 0 Such pathological tissue 
IS also found to contain more water It would thus seem from 
the tissue characteristics that an imaging method based on 
water content and/or would easily differentiate between the 
various tissue types yielding cross-sectional images showing 
anatomical detail as reflected in the distribution of vjater 
throughout the body The spin-spin relaxation time T 29 which 
also varies with tissue type, might be similarly used, given 
suitable modification of the imaging method Imaging capabi- 
lities of these two important parameters, and T 29 together 
with the spin densities make NMR imaging a unique, versatile, 
and powerful technique in medical imaging 

1 2 BASIC NMR IMAGING COMPONENTS [ll] 

A typical NMP imaging system designed for human imaging 
IS shown in Fig 1 4 [7] The basic comoonents of an NMR 
imaging are as follows 

1) A magnet, which provides a strong (range - 500 G to 15 KG), 
uniform, static magnetic field (H^) 

2) An RF transmitter, which delivers radio frequency (range - 
2 to 70 MHz for proton imaging) magnetic fields to the 
sample (Hj^(t)) 

3) A gradient system, which produces time varying magnetic 
fields (range - 0 01 to 1 G/cm) of controlled spatial 
nonuniformity (G^(t), Gy(t), (j^Ct)) 


8 
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4) A detection system, which yields the output signal 
(y(t)) 

5) An imager system, which reconstructs and displays the 
images 

Imaging system shown in Fig 1*4 uses a single RF coil 
which IS used for both transmission of RF signal (longest RF 
pulse duration is on the order of 2 ms) and reception of FID 
signal (like a Pulsed Radar) The stochastic imaging dis- 
cussed in the next section assumes two coils one for trans- 
mission and other for reception, since transmission and 
reception are done simultaneously 

Resolution in a conventional NMR imaging system depends 
on tho field homogeneity and available strength of field 
gradients For example, with a magnet of the static magnetic 
field of = 5 KG having homogeneity of =4 ppm 
(0 02 G) with field gradient of ^ = 04 G/cm, the 

expected resolution limit is approximately 0 5 mm The physi- 
cal arrangement of components and instrumentation is an 
important factor in NMR tomography 

1 3 HISTORICAL BACKGROUND [l8] 

At the beginning of 1970’ s both Lauterbur [l6] and 
Damadian [8] showed that NMR spectroscopic techniques can be 
applied to imaging of the human foody and that the techniques 
have diagnostic uses similar to those of X-ray computerized 


10 


tomography (CT) In 1978, Andrew et al demonstrated the 
very-high-resolution capability of NMR by obtaining a fine 
detail imago of the submillimeter septum of a small orange 
More rocently, human head images shown by Moore and Holland 
et al have definitely demonstrated the great potential of 
NMR imaging in medical imaging Already the subject has 
acquired a variety of names of which more common are NMR 
imaging, spin mapping, Zeugmatography , NMR tomography 

1 4 SYSTEM THEORY APPROACH 

It is well known that a linear time invariant (LTI) system 
IS completely characterized by its impulse response System 
identification problem is equivalent to the impulse response 
determination Various methods used in impulse, response 
determination of LTI systems can bo broadly classified as 
sinusoidal wave testing, pulse testing and cross correlation 
method These are the familiar techniques already in usage in 
NMR In correlation method an input signal with impulsive 
autocorrelation function is applied to the system and the 
response is crosscorrelated with shifted versions of input 
signal to obtain impulse response point by point 

In case of linear timo varying (LTV) systems, response of 
tho system for an arbitrary input cannot be given in a closed 
form as in the case of LTI systems But in general when infi-* 
nite series solution is obtained under conditions which permit 



11 


the system to be viewed as under perturbation, the solution 
resembles that of a LTI system In NMR imaging equations 
governina the system behaviour (Bloch equation) correspond to 
that of LTV system under the conditions, the input RF excita- 
tion IS small [12] the response resembles the familiar convo- 
lution integral and hence emenable to cross correlation 
approach The output of cross correlation operation gives 
the information regarding the spatial spin distribution [2] 
Details of this approach are given in Chapter 2 

1 b SCOPE OF THIS THESIS 

Chapter 2 develops the mathematical model of NMR imaging 
problem, NMR imaging is treated in the system theoretic point 
of view as LTV system, described by state representation Solu- 
tion of the state equations is discussed Chapter 3 covers 
the cross correlation techniques to process the FID signal to 
obtain the spatial spin distribution Chapter 4 discusses the 
computer simulation of NMR imaging system model and the signal 
processing schemes FID signal corresponding to specified 
gradients and RF excitation are generated by computer using 
the model Associated cross correlation scheme is applied to 
obtain the assumed spin distribution Chapter 5 concludes the 
thesis with suggestions for future work 
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CHAPTER 2 

MATHEMATICAL MODEL 


2 1 BLOCH EQUATION 

The Bloch equation [l], gives an accurate but phenomenolo- 
gical description of the time dependence of the ’nuclear magne- 
tization’ M(t) in the presence of an applied magnetic field 
H(t) The nuclear magnetization M(t) is the source of the *NMR 
Signal’ from which the image is ultimately constructed The 
NMR sample i$t re ate d as a black box with H(t) as the input signal ? 
or stimulus, and M(t) as the output signal or response The 
black box is characterized byM^Cxjyjz), and T 2 ? and its 
behaviour by the Bloch equation 


dM 

dt 


_ _ M u + M M - M 

= YM X H - -T - -y-y -) - (-A= -) u 


(2 1 ) 


2 n 

The above equation can be conveniently represented in matrix 
notation as 










1 

yh 0 



“y 


- yh 

- ^ YHi(t) 



L ^ 


0 

- yHj^(t) 0 


L 2 j 
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List of Symbols 


f 

Im 

Re 

H 

o 

Y 

h(t,x,y,z) 


H(t,x,y,z) 




M(t,x,y,z) 


M = M^+3M 

“o = “YHq 
10 

U^.Uy.Uz 

G^(t),G (t), 
G,(t) 


Three dimension Integral 
Imagionary part of 
Real part of 

Strength of the static magnetic field 
Gyromagnetic ratio For protons its value is 
4 26x10"^ Hz Tesla"^ 

Gradient magnetic field 
(G^(t)x -f- Gy(t)y + G^(t)z) 

Its effect in Uy,u„ is neglected because u term 
dominate s 

Total magnetic field experienced at t|X,y,z 
(represented as H or H(t)) 

It has components in Uy? u , , and u„ directions 
Input RF signal 

Its effect IS only Uy direction 
Sample magnetization at t>x,y,z 
(represented as M or M(t)) 

Components of M 

(Complex) transverse magnetization 

Equilibrium magnetization or spin density 

Larmour frequency corresponding to the static field 

Irradiation frequency of the RF field 

Unit vectors in x,y,z directions respectively 

Gradient magnetic fields applied along x,y,z 
axes respectively 
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The nuclear magnetization M(t) , induced in the sample by 
the magnetic field H(t), is the local sum of the magnetic 
fields of the protons The magnetization is a bulk property 
of the sample rather than that of individual protons 

2 2 EXCITATION METHODOLOGIES 


Two kinds of excitations are seen in NMR imaging First j 
short-pulse excitation in which first a RF pulse, sufficient 
to tilt magnetization into transverse plane, is applied At 
the end of RF pulse, field gradients are applied during which 
FID signal is recorded Second, both excitation and gradients 
are applied simultaneously throughout the experiment and FID 
signal is recorded 


Short Pulse excitation method can be conveniently explained 
through the Bloch equation discussed earlier The system 


dynamics after the application 

of excitation pulse is 
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to 


r 


(x,y,z) 


M (x,y,z) 
^0 


0 


— ^ 

where M^lx^y^z) = V (x,y,z) + M (x,y,z) 

0 ^ o 

M can be given as 


Solution for 


t 

-t/T2 -JY / h(t')dt' 

M = e e o M (x,y,z) 

where M (x,y, 2 ) = M (x,y,z) = -jM (x,y,z) 

It IS interesting to note that similar solution is arrived at 
in stochastic NMR imaging 


2 3 PERTURBATION ANALYSIS [12] 


In this method both input excitation and gradients are 
considered at same time 

Equation (2 1) can be written as 


dM, 

dT 


M. 


= - ^ + M vh 


6U 

dt " 


M. 




dM, M, - 

dT = - \ - (-S7- ) 


Let M = + jM 


(2 2 ) 

(2 3) 

(2 4) 
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The equations (2 2) , (2 3) and (2 4) may be written, neglecting 
spin-lattice relaxation (i e , t << T^) in the form 

If + (f~ + JYh)M = JYH3^(t) (2 5) 

dM„ 

fy! -yH^Ct) Im M (2 6) 

Due to linear time varying characteristics of a driven NMR 
system it is usual to discard the simple methods of Fourier 
analysis applicable to other spectroscopic forms The parti- 
cular case considered here is the response of the NMR system 
to an arbitrary shaped pulse of duration less than T 2 The 
approach to be employed is a combination of perturbation theory 
and linear system analysis [19] 

Let us recall that M is a complex quantity, M +jM , and 

X y 

real These can be expanded as 

M = f + (2 7) 

and are small quantities of the nth order Assuming 

that the driving function Hj^(t) represents but a small perturba- 
tion to the Hamiltonian of the NMR system, 1 e , applied exci- 
tation amplit-ude is far less than the static magnetic field, and 
that and are zero when Hj^(t) =0 Substituting eq. 
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(2 7) in (2 5) and (2 8) in (2 6) yields a set of equations 
of various orders, and representative equations of one order 
are 




+ ( 


dt 




oi 


z 

dt 

asL 

Thus , 


may 

terms 

of M 

n-2) 


(n) ^ 


(n-1) 

z 


(n-1) 

‘z 


, M 


(2 9) 

(2 10 ) 

(n-l) 


and so on In other words, knowing the 
lowest-order terms, all other terms may, in principle be 
calculated from these recurrence rela tionships Quan titles 
such as YH^(t) are of order n as is small 

Consider first the zeroth-order terms 


dt 


+ ( 4 - + jYh) = 0 
*2 


(2 11 ) 


However, in the absence of any irradiation, we know that M 
must be zero and so =0 (No tilting of magnetization 

vector into transverse plane and hence no reading in that 
plane) The equation in M gives 


dM 

Ht 




(2 12 ) 


and so is constant which may equate with the equilibrium 

magnetization M^(x,y,z) Moving on to the first-order equa- 


tions 
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dM 

Ht 


( 1 ) 


+ (~ + jYh) = 3YH,(t) = 3YHi(t) M^{x,y,z) 


o 

(2 13) 

This equation will be analyzed in considerable detail later 
Also, 


dM 


(1) 


r\} 


- — = 0 


( 1 ) 


This IS a repeat of eq (2 12), but as ' is zero for Hj^Ct) = 0, 
it must be zero all times 

( M = I- , 

z 

~ +• O 4* ) 

The second-order equations thus give 
dM^^^ 


dt 


4- + jyh) = 0 ( = 0) 

* rs Z 


( 2 ) 


and here M' ' is zero for the same reasons 


dM 


(2) -L 


dt 




Hence 

t 

tSS f -YH,(t) dt 

•-O0 y 

This term describes the departure of the magnetization from the 
equilibrium alignment along the z-axis Finally, the third 
order equation in M is 
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+ (™ + jYh) = jYHi(t) 


and this expression describes the deviation of M from the 
linear behaviour of q (2 13) 

Equation (2 13) describes a linear differential equation ^ 
and so the technique of linear systems analysis may be brought 
to bear on the problem of solving tho equation [19] Hence, 

- / (^ + 0Yh)dt' ^ / ^(y +JYh)dt' 

mO) = e ° 2 j MgCx.y.z) e° ^ dtj_ 

o 


(2 14) 


writing Eq (2 14) in detail we have 


M 


■j. 

=r oxp [-/ ~ + JYHQ+JYxG^(t‘)+JYyGy(t*)+ 0 YzG 2 .(t*) dt’] 


0 "2 
t 


/ JYHj^(tj^) M|j(x,y,z) exp [ f y“ ■^JYHo+3YxGx( t' ) 

o 2 


+ OYyGy(t') + jYzGz(t') dt'] dtj_ 


(2 15) 


Substitution of =5X yields 

t 

e 
o 


M 


(1) = -jYMQ(x,y,z) / Hj_(t~T) e*" 

2 


t-t: 


t-x 


t-T 


JYX /G^(t')dt' „JYy / G„(t')dt' „JYZ / G,( t' )dt' 


t 


(2 16) 
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With (2 lo) as the base, soiuti on for rerndining terms can 

be written, at least in principle Hoult [l2] showed in 

computer simulation studios that under the assumption that when 

input excitation is small enough to tilt magneti .sation vector 

from static field alignment by an angle less than 30° the 

( 3 ) 

effect of and higher order terms can be neglected and the 

dominant term is only Our further analysis takes only 

( 1 ) 

term into consideration 

Hence, receiver coil output is spatial integral of 
term only and estimation of spatial inf ormation( spin density) 
on this output IS a fair approximation 

Similar solution (Eq (2 16)) can be arrived at using 
Kinariwala’s approach [lO] This approach is to decompose the 
square matrix represented under Eq (21) into two matrices, 
A^{t) and Aj^(t)j whore AQ(t) satisfies the commutative condi- 
tion A 2 ^(t) is interpreted as a perturbation upon AQ(t) This 
method is very similar to Hoult* s perturbation analysis, except 
whole matrix is considered instead of individual terms 
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CHAPTER 3 

signal processing 

3 1 NMR DETECTION SYSIBA 

Expression for M(t), the nuclear magnetization is 
derived in the previous chapter It is the function of the NMR 
detection system^ the receiver, to detect iS(t) and generate 
an output signal y(t) The receiver coil, which usually 
surrounds the samole, is an antenna which picks up the fluctua- 
ting nuclear magnetization of the sample and converts it to a 
fluctuating output voltage V(t)? where 

V(t) = - ^ / M(t,x) B^(x) dx 

The function B^(x) describes the sensitivity of the receiver 
coil at different points in space The primary objective of 
receiver coil design is to prescribe wire placements so that 
B^(x) has the largest possible transverse cornoonent The 
longitudinal cornoonent of B^{*x) contributes little to the out- 
put voltage, and can be ignored This is a result of the fact 

that the time derivative of is much less than that of 

transverse component M 2 ,(t,x) decays exponentially with the 
time constant T^ typically 01 to 1 S, while the transverse 
component is oscillatory with a period of typically, 0 Ob to 
0 2 ps Let 
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B^(x) - au + bu„ 

where a and b are fixed, but unknown constants 

Instead of V(t), time integral of V(t) (let this be y*(t)) 

IS considered which leads to mathematical simplicity Hence 

y'(t) =7 (aM^(t,x)4- b My(t,x)) dx 

= f Re((a'-jb) M(t,x)) dx 

Making the replacements a - I<cos0' and b = -Ksin0' 
y’ (t) ~ / Re ICe M(t,x) dx 

= - / Re K M(t,x) dx 


where K is complex arbitrary constant 


Dropping the constant part we have f M(t,x) dx, and let this be 
represented by y’*(t) In detail y”{t) can be written from 
Eq (2 16) as 


y” (t) 


t ^ jyH ) 

= fff M (x,y,z) / H, ( e 2 

o 

t— t— ^ 

/ Gjj(t')c3t' gJYy / G (t')dt' 
t t ^ 


t-r 

gJYz / G^(t')cit' 


By quadrature phase detection term can be eliminated 

Let y*' (t) after quadrature phase detection be represented as 
y(t) 7^1 th the definitions 
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t-T 

u(t,'’^ ) = Y / ') 

t ^ 

v(t,T) = Y / Gjt') dt' 
t ^ 

and 

w(t,T ) = Y / G^(t') dt' 

t 

y(t) can be written as 

+ «> t -Vt^ Jxult^T) 

y(t) = fff Mq{x,y,z) / H^(t~t ) e ^ e 

— oo O 

e3yv(t,r) gjzw(t,r) ^3 

The imaging experiment is treated in terms of system 
analysis The underlying concept is that of a black box which 
IS tested with an input signal H^(t) and responds with an out- 
put signal y(t) The feature to be revealed in this imaging 
experiment is the spin density distribution of the biological 
sample as a function of space coordinates For and T2 
imaging a different mothodology has to be adopted From Eg 
(3 1) it can be seen that four variables (shape of the RF input 
signal, and the modulations of magnetic x,y,z gradients) are at 
our disposal Consequently the imaging experiment may be treat' 
ed as a four-input timo invariant system or a one-input, time 
deoendent system, when tho time dependence is imposed by the 
magnetic field gradients and the one input is the RF signal 
More specifically, the NMR imaging experiment is described in 
terms of a linear time-dependent system 



Broadly imaging modalities may be classified as being 
'reconstructive' or ' non-reconstructive’ The former technique 
IS an adaptation of the principles of image generation in X-ray 
CT, and involves reconstruction of the image in its entirety 
from integral progoctions 

Modalities in which an imago reconstruction scheme is not 
involved are classified as non-reconstructivc These can be 
classified in two different ways 

(i) Based on dimension 1-Dimcnsion, 2-Dimension, 3-Dimension 
Using appropriate excitation signal and gradients, it is 
possible to excite desired section in spin density 
volume (Fig 45) 

(ii) Based on pulse sequence Van^ies of pulse sequences are 
in usage A specific method known as Kumar, Welti, and 
Ernst (KVVE) direct Fourier NMR tomographic imaging is 
discussed below 

The RF pulse and gradient sequence used in KWE method is 

illustrated in Fig* 3 1 and commences at t = 0 with a short 

nonselectivc RF pulse to nutate some or all of the equilibrium 

magnetization into the x,y plane This is then followed by 

periods t and t of variable length during which linear mag- 
X y 

netic field gradients 0 ,G are applied in x and y directions 

X y 

respectively Finally, the FID is sampled in the presence of a 
z-gradient for a time t^ = 
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ThG z component of the magnetic field is therefore given 
by 


Ho + G^x 




“o ^ 7 ?- 


0 < t < 

< t < t^+ty 


and the response of the spin system y(t) will be a function of 
and t , as well as i e , y(t) = y(t^>t^,>t^) The cxperi- 
ment is repeated for a full set of regularly stepped or mere 
mented t^ and t^ values allowing a three-dimensional image of 
the object to be reconstruc ted as discussed below 


The observed signal y(t) can bo written from Eq (3 1) as 

_ ... JYG^xt JYG yt oyG^zt -( 

y(t) = JJJ M^(x,y,z) e 0 ^ e ® 

dx dy dz 

Neglecting the effect of T2, i\/tQ(x,y,z) can be obtained by 

inverse 3-D Fourier transform as given below 

... -OYG x«t -jyG^y't 

M^CxSy^z') = /// ^ e ^ ^ e 


di^ dty dt^ 


3*2 STOCHASTIC NMR IMAGING [2] 


Stochastic excitation denotes a noise-modulatod RF input 
signal, which is supplied for the entire length of the imaging 
experiment RF noise is conveniently generated by a sequence of 
RF pulses with random flip angles The pulses are separated by 
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the sampling interval At For each input pulse one value of 
the response is sampled The stochastic RF excitation delivers 
its excitation energy to the NMR system within a characteristic 
relaxarion time While pulsed excitation concentrates the 

same energy into a single pulse of length ^t 

The spin density is retrieved from cross correlation of 
input and output records The longer these records, the better 
will be achievable signal to noise ratio (SNR) of the spin 
density In stochastic system analysis the input signal usually 
IS white noise with zero mean so that the autocorrelation of the 
real-valued input is proportional to a delta function 

T 

Lim i S H, (t) Hn(t-a) dt = (o) (3 2) 

T-^ ^ o 

The proportionality constant ^2 second moment of the 

power of the input 

For time-invariant systems the linear impulse response 
function can be retrieved from cross-correlation of input and 
output signals according to 

1 T 

Lim Y S y(t) Hj^(t-o) dt = ^2 i^) 3) 

T-> o 

where M‘(t) is impulse response, 

_t 

M‘(t) = S M^(x) exp [- |- + 3YX / Gy(t’) dt‘] dx , 

^ ^2 o ^ 
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with M'(t < 0) = O The cross-correlation is a time average 
of the product of the system response with the delayed input 
signal The result is the linear impulse response as a fun- 
ction of the delay o In case of time-dependent systems this 
concept cannot bo applied However, for the imaging system 
the goal is to retrieve its time-invariant part but not the 
complete time dependent imoulse response Since the time 
dependent part is known, the system response y(t) can be 
correlated with the product of ail the time-dependent factors 


in the response integral 

-jx' u(t,cf) -oy’ v(t,-cj) w(t,e) 

P(t-o) = e e e Hj_(t-o) 


(3 4) 


The primed coordinates as well as o become variables in the 
cross-correlation Its evaluation yields the spin density as a 
function of the cross-correlation variables* The complete 
scheme is shown in Fig 3 2 [3] 

Lim 4 / Y(t) P(t-o) dt 

oo * o 


oo -r/r T jxu(t,T: ) 

= SSS Mq(x,y,z) / e ^ Lim j f ^ 

0 T oo o 

-jx'u(t,c?) jYv(t,T) -jy'v(t,o) -jz'w(t,a) 
e e o e 

gjzw(t,t:) 


Hj^(t-T) H 2 _(t-a) dt dr dx dy dz 


(3.5) 
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= (2it)^ (12 /// MqCxjYjz) / e 


-■C/T2 


»(x-x') ft(Y-y') ■tj(z-z') 


•15 ( r -Of) dT dx dy dz 


(3 6) 


= (2-jt)^ (i 2 Mq(x' , y' »z' ) 


(3.7) 


In Eq (3 5) the cross-correlation has been written down 
explicitly but the order of integration has been changed The 
stop from (3 5) to (3 6) involves the evaluation of the inner 
most integral The fourfold product of delta functions 
result from Eq (3 2) , the definition of the delta function 


1 

Sc 


CO gk(x-x 
/ G 


= «>(x-X^) 


(3 8) 


and a postulated orthogonality of the four functions in the 
product of Eq (3 4) » The step from Eqs {3 6) to (3 7) is 
the evaluation of the delta functions 

For the orthogonality to be satisfied three conditions 
must be met 

Condition 1 

The input signal Hj^(t) must have zero mean 

T 

(1, = Lim 4 / H, (t) dt = 0 (3 9) 

oo o 
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Condition 2 

The time integrals over the magnetic field gradients 
must always bo different from zero 

u(t,T), v(t,x), w(t,r) 5^ 0 (3 10) 

for all t and some given cr 
C ondition 3 

The difference of the local phase angles in the different 
directions x,yyZ should never vanish for all primed and un~ 
primed values of x,y and z 

xu(t,a) - yv(t,a) 0 , 

yv(t,a) - zw(t,a) ^ 0 j (3 11) 

zw(t,o) - xu(t,a) 0 y 
for all times t, the o' given, and all x,y,z 

Condition 1 is trivial Condition 2 can be fulfilled by 
modulating a field gradient such that no gradient component 
ever changes sign The last condition is the most restricting 
one One way to observe it is by using significantly different 
field gradient in each direction x,y and z A different gra- 
dient scheme is simulated in noxt chapter for the case where the 
spin density shall be recovered at discrete coordinates y‘ , 
and z’ 
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The final result of the cross correlation (3 7 ) is the 
spin density as a function of space coordinates, scaled by 
the input power, p.2# 3 negative exponential involving the 

transverse relaxation time T2 and the cross-correlation delay cr^ 
Input power 1I2 cannot be made arbitrarily large in order to 
increase the system response, since eventually the system will 
become nonlinear and the response will be damped However, 
there exists an input power levol which yields a response maxi- 
mum and consequently optimum SNR in the acquisition of the res- 
ponse 

In comparison with pulsed FT NMR, stochastic resonance [4] 

1) has the same sensitivity , 

2) allows an independent optimization of sensitivity and 
resolution in contrast to pulse FT NMR $ 

3 ) reduces the problems in the receiving electronics associa- 
ted with the large dynamic range of impulse signals, 

4 ) reduces the peak excitation power such that wider dynamic 
ranges may be covered , 

5 ) admits coherent averaging of the system response in conne- 
ction with cyclic excitation , 

6 ) will exhibit different saturation features In particular, 
the nonlinear response will involve double resonance effects 1 

7 ) investigates the spin system in dynamic equilibrium with 
the excitation and not the unperturbed thermodynamic equi- 
librium state 
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Stochastic excitation in NMR imaging has three advantages 
over conventional methods with sparsely pulse PF input 

1) The excitation power is reduced by several orders of magni- 
tude 

2) Since no deterministic pulse cycle exists, every addi- 
tionally acquired data point increases the quality of the 
entire image 

3) For slowly moving objects, reconstruction of the image can 
be restricted to those data, during the acquisition of 
which the object was in approximately the same state 

33 CHIRP IMAGING 

In NMH imaging the quest for high spatial resolution is 
of fundamental concern in clinical applications High resolu- 
tion necessiates the use of strong magnetic field gradients 
Conventional imaging schemes traditionally employ fast switch- 
ing of gradients Since the spt,ed of switching imposes an 
upper limit on the usable gradient strength, it is desirable 
to search for imaging schemes which permit slow switching of 
field gradients admitting stronger gradients for improved 
resolution NMR imaging with sinusoidal field gradient modula- 
tion [5], [9] seems to be an attractive solution# 

In general an infinite number of gradient modulation fun- 
ctions (G^(t,T), G,,(t,X)> G_( t,^ )) are conceivable, which 
obey the three conditions, Eqs (3 9), (3 10) and (3 11) We 
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propose a new class of imaging scheme, viz , Chirp processing 
because of its simplicity in generating gradient functions and 
subsequent processing structure Chirp processing has been 
dovoloped to provide a solution for tho conflicting require- 
ments of simultaneous long-range and high -resolution perfor- 
mance in radar systems [l4] Chirp processing is often inter- 
preted as an alternative to Fourier transform in signal proce- 
ssing 

Consider Eq (3 l) After the application of ix/2 pulse, 
y(t) can be written as 

“t/T 2 3 xu( t) jyv( t) jzw( t) 

y(t) = fjf MQ(x,y,z) e e o e dx dy dz 

(3 12) 

T 2 IS a spatial function and hence cannot be brought outside 
the integral, but let us neglect the effect of T 2 

Let the gradients bo linear functions of time as given 

below 


G^(t) 

= K^t 

(3.13) 

QyCt) 

= K2t 

(3 14) 

Gj,(t) 

= Kgt 

(3 15) 


2 

t K, t^ 

u(t) = / K, t' dt' = — 

0 


and 


(3 16) 
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Similarly v(t) = (3 17) 

2 

w(t) = — (3 18) 

where ^3 constants The three conditions Eqs 

(3 9), (3 10) and (3 11) are satisfied by this particular 
choice as per the constraints given below 

1) Input signal H^(t) must have zero mean 

Tp 

/ Sin o)t dt = 0 
0 

2) The time integrals over the magnetic field gradients must 
always be different from zero 

Our choice of field gradients ^which are linear fun- 
ctions of time satisfy this property (Eg (3 13), (3 14) 
and (3 15)) 

3) The difference of the local phase angles in different 
directions x,y,z should never vanish for all primed and 
unprimed values of x,y and z From Eq (3 11) 

X u( t) - y v(t) 0 

X ^ y ^ 

This IS possible +o maintain except for some set of values 
of x,Kj^,y,K 2 All points in x-y plane can be imaged except 
those lying on straight lino y = This problem can be 



35 


overcome by repeating the experiment with different Chirp 
rates, i c , and instead of and K 2 such that 
Ki/K^ ^ K1/K2 

The constraints under this condition can be given as 


i yKg -4^ zK^ 
and x’ k, 4 y'k^ 4 z' 

-3X'K,t^/2 -jy'K2t^/2 

Multiplying y*(t) (Eg (3 12)) with e o 

-JZ'Kgt^/2 

e and integrating with respect to time amounts to 

jK2(y-y')t^/2 J'<3(z-z')t2/2 

/// MQ(x,y,z) / e e 

dt dx dy dz 

In the inner integral each term behaves approximately like 
a delta function 

00 jK. ( x-x' ) t^/2 

Consider d(x-x') = / o dt This can be 


written as 


t^/ 20 ^ 


where o - 


/ e 
1 1 


\/k, (x-x* ) e 


31^/4 


From the knowledge of Gaussian function this integral can be 


written as 
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, neglecting the phase part 


The above function behaves approximately like a delta function 
and be represented as d(x-x') When x = x’ the denominator goes 
to zero and hence the total quantity to When x* ^ x then 
the function value is finite, ideally speaking it should be 
zero x’ >> X of X* << x the function becomes zero But such 
cases are not encountered due to finite dimensions of the 
specimen 

Since conditions 1,2 and 3 are satisfied the four inte- 
gral function reduces to K" M^(x‘ ,y* ) , where K” is a con- 
stant Vi/ith different x*,y*,z^ the whole specimen can be 
imaged Simulation work carried out in Chapter 4 justifies 
this method The processing can be made real time, at least 
in theory with an array of chirp multipliers and integrators 
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CHAPTER 4 

SIMULATION 

4 1 STOCHASTIC NMR IMAGING SIMULATION 

The excitation Hj^(t) can be realized with either a 
white sequence of pulses or continuous noise with zero mean 
The particular distribution function is of no influence for 
evaluation of a linearly driven system In particular a 
binary sequence of flip angles will be sufficient Such 
sequences are easy to generate with feedback shift registers 

[15] 

4 11 Discretization of the Spin Density 

The second orthogonality condition is fulfilled by 
choosing distribution functions for the field gradients 
G^>Gy,G^ with non zero mean 

In consideration of the third orthogonality condition 
it will be assumed for now that the measured response derives 
from a spin density, which, contrary to reality, is defined 
only on a discrete set of equidistant coordinate oointsx,y,z 
This set of points shall be the same as the one on which the 
spin density Mq(x' , y’ fZ* ) will be reconstructed by cross- 


correlation 
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To avoid overlap of values, an interleaving set of 
discrete coordinates is constructed by setting 

X' = (n^+ oc^) a , 

y' = (ny+ tty) a , 

z' = (n^^ tty) a , 

where n^ is an integer, 0 < < 1 for i = x,y,z and the 

dimension of a discrete spin density element or picture cell is 
apecifiod by 'a* Similarly a discrete range of values is 
introduced for the gradient integrals 

U = 

V = (n^ i- ^pg 

where is an integer, 0 < < 1 for i = u,v,w and the 

product of the minimum field gradient step and the gradient 
switching interval ^t is denoted by g Without loss of 
generality the field gradient switching interval will be 
assumed to be identical to the sampling interval, and hence 
the same symbol The unit of g is rad/m The value of g 
measures the minimum precession phase per picture ceil which 
IS acquired during one gradient step of length At For a 
symmetrical distribution of gradient numbers n^, the mean 
field gradients are g^ = <x ^ g/ At 
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Evaluation of one. of equations (3 11) provides insight into a 
choice for cc ^ and 

x'u - y'v = [(n^+ ax^^’^u'*' “u^ " “y)("v’^ ® 9 

The orthogonality condition states > that the right-hand side 
must bo different from zero for all values of 
n^^ny^Hw Setting = a^j = ay^ it can be suen 

that 

x'u - y'v = l^'^x'^u'*'^"x'^*^u^ “x"^® ®V ““y^® ^ 

With 

“x “u ^ 1/7 ? 

ay - ay = 3/7, {4 1) 

the third orthogonality condition is fulfilled by the inter- 
leaving coordinate and gradient integral values according to 
Cqs (3 10) , (3 11) , and (4 l) 

412 Numerical Simulation 

The linear response has been calculated according to 
Eqs (3 1), (3 10), and (4 l) The spin density was defined 
on a 3x3x3 mesh (Fig 4 1) with values reflecting proton 
densities of living tissues The transverse relaxation time 
was T^ = 0 1 sec (As mentioned earlier this should be a 
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spatial function, but hero it is assumed to be constant) To 
keep the computation times at an acceptable level, the product 
of spatial resolution and precession was set to ag = 2 b rad 
Tho sampling interval was 0 01 sec For each sampling interval 
new field gradient values and a new excitation value were 
derived from a 20 stage binary feedback shift register The 
excitation could assume two values, = [~l/2, +1/2] with 
equal probability, while the field gradient numbers could 
assume four values, n„, n.., n^ = [-3/2, -l/2, +i/2, +3/2] with 
probabilities l/S, 3/8, 3/8, and l/8 The field gradient num- 
bers were obtained by shifting the feedback register thrice and 
adding the binary outputs of each shift 

The cross-correlation was executed for a 5x5x5 point mesh 
of the spin density M (x',y*,z’) and a time shift o = The 
result is shown in Fig 4 2 in terms of slices showing y‘ z’ 
planes for different values of x’ The numbers denote the 
retrieved spin density values percent The background noise 
results from the finite sample record length 

For a random data (0 0 < M^(x,y,z) < 1 0) also the simu- 
lation is earned out Tho assumed spin density and imaged 
values arc shown in Fig 4 3 All tho spin densities shown in 
this chapter are percentile 
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Fig 4 2 Slices of retrieved spin density 
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Fig 4 3(a) Spin density used in random simulation It 

IS defined as a 3x3x3 array of cubic picture 
cells with dimension a 
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Fig, 4 3(b) Retrieved Spin Density Slices 
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4 2 CHIRP IMAGING SIMULATION 

Under this imaging scheme, wo nssume that 7i/2 pulse is 
applied and then FID signal is recorded in the presence of 
linear time function gradients The pulse sequence is shown 
in Fig 4 4* 


4 2 1 Simulation 

volume 

This process assumes that only a line of wholo^is excited 
and FID signal from other part is zero (Fig ^ 5(b)) For one- 
dimension chirp imaging the FID signal can be written from 
Eq (3 12) as 


y(t) 


r , s F /2 

= / M^(x) e dx 


Here T 2 offoct is not taken into consideration Vihen T 2 is 

assumed to bo a constant its effect can be taken outside by 

multiplying y(t) with 65536 samples of y(t) are multi- 

-j X* Kt t^/2 

plied with o and av^rigedover time to get M^(x‘) as 

shown below 


/ y(t) G 


-jx'K3_F/2 


dt 


jK, (x-x') t^/2 

= / Mq(x) f 0 ^ dt dx 


= f Mq(x) d(x-x' ) dx 

= Mj^(x') 



ke 


90 ° Pulse 


90° Pulse 



Fig 4*4 Pulse sequence for 2~din'ensional chirp 
imaging 


S3. 3^ ^ 3 



Fig 4*5 


cheme of classification of Imaging techniques ( a) sequential 
oint measurement (b) sequential line measurement , C^) sequential 
lane measurement, and (d) simultaneous volume measurement 
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For simulntion, Mq(x) is takon as 


Mq(1) = 19/ y M^(2) = 73/, M^(3) = 60/, and imaged values are 
(using piograms 5 and 6) 

Mq(1) - 25/, Mq( 2) = 77/, M^(3) = 64/ and h\^{ l) == 5 6/, 

Mq(4) = 6 6/ 

Ideal values of and M^(4) should bo zero Those two 

values givo idea of noiso? where spin density is not present 


422 2-D Simulation 


This process assumes that only a slice of whole volume is 
oxcibed [6], [20] and TID signal from remaining part is zero 
(Fig 4 5(c)) For two^ dime ns ion chirp imaging the FID signal 
can be written from Eq (3 12) as 

jK,xt^/2 jlC^xtVs 

y(t) = // MgCx.y) e o dx dy 


65536 samples of y(t) are multiplied with o 


-jKnX' t^/2 


^ and averaged over time to get M(x' given 


o 

below 


-jK3_x-t2/2 -3l<2Y't^/%^ 


/ y(t) 0 

jK,(x-x' )t^/2 . A 

•. V Q ^ dt dx dy 

= // M„(x,y) Jo o 


= // M^{x,y) d(x-x') d(y-y') dx dy 
= Mq(x' ,y' ) 
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^^(xjy) IS tnken as slices (Fig 4 6) of spin density shown in 
Fig 4 1 The imaged values (Programs 7 and 8) are shown in 
ng 4 7 

Programs 7 and 8 are ran for random 2-D (slice) spin den- 
sity Results are given in Fig 4 8 The first lino in this 
slice IS the same data ran for 1-D chirp imaging (Section 
4 2 1) The iinaged values for the first line are different 
from 1-D chirp imaged values which clearly indicates the 
correlation effect of two impulses, viz , a(x-x') and d(y-y’) 

4 3 Dc. TAILS OF PROGRAMS 


The simulation work is earned out on DEC-1090 computer 
system Storage of 65536 comolex values of FID signal re- 
quires about 2800 blocks The following gives details of 
programs The listings are given in Appendix 


Program 1 

Program 2 

Program 3 
Program 4 
Program 5 
Program 6 
Program 7 
Program 8 


Simulates stochastic magnetization - 
CPU time 6 hrs 5 mts 24 55 secs 
Cross-correlation for stochastic mR tomography 
0 34 50 72 

Formated printing of input 
Formated printing of output 
Simulates 1-D chirp magnetization 

1- D chirp imaging 

Simula tGs2-D magnetization - 0 2 53 12 

2- D chirp imaging -07 9 36 
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CHAPTER 5 

CONCLUSIONS 

NMR techniques can be applied for imaging the protonic 
tissue matter in the whole body The image can be obtained 
by a suitable set of gradient and RF signals and performing 
the appropriate signal processing The image can be obtained 
point by point or line by line or slice by slice or total 
volume Notable among the schemes are Fourier transform (FT) 

techniques and the correlatioh techniques We have used PN 
sequence as the RF signal and spatially linear gradients as 
suggested by Blumich The FID signal was simulated on the 
computer for a 3x3 matrix of 3 slices of spin density The 
simulation of the FID signal took about 6 hours CPU time on 
DnC-1090 computer, and the processing on FID signal to recover 
tho image about 30 minutes Due to this large computer times 
involved larger matrix sizes, larger number of slices could not 
be tried out A new schomo of obtaining the image through the 
chirp processing is implemented In this scheme an RF signal 
is applied for a short duration followed by application of 
linear spatiaJly and temporally varying gradient fields The 
FID signal starting from the point of termination of RF signal 
is simulated, and tho associated chirp processing is carried out 
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ThG reconstructed images were found to be in good agreement 
with original spin density The chirp method is applied for 
1~D and 2-D imaging satisfactorily But the results for 3-D 
wore not encouraging 

For obtaining the NMR image of 3x3x3 points from its FID 
signal aoproximately 1,769,472 complex multiplications and 
additions havo to bo performed by computer, which require a 
largo computational time The time requirea under the modali- 
ties involving FT techniques can be significantly cut down 
using Fast Fourior transform algorithms HovJever, chirp 
processing cannot be made faster quo to lack of fast algori- 
thms It appears possible in any case to reduce tho storage 
roquiroment and CPU time by using sories-parailel processing, 
and should be attempted 

A detailed examination of the 3~D chirp is required to 
identify bho reasons for unsatisfactory results we obtained 
We havo not apolied any windows on tho data to study their 
effects on side lobe levels It appears possible to use 
appropriate window to improve the side lobe characteristics as 
done in conventional signal processing on finite data 
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